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Abstract 

We investigate in the Hartree Fock approximation the temperature 
and filling dependence of the superconducting 7r-phase for the Penson- 
Kolb-Hubbard model. Due to the presence of the pair-hopping term, 
the phase survives for repulsive values of the on-site Coulomb interac- 
tion, exhibiting an interesting filling and temperature dependence. The 
structure of the self-consistent equations peculiar to the 7r-phase of the 
model allows to explicitly solve them for the chemical potential. The 
phase diagrams are shown and discussed in dimension 2 and 3. We also 
show that, when a next-nearest neighbours hopping term is included, 
the critical temperature of the superconducting region increases, and 
the corresponding range of filling values is shifted away from half-filling. 
Comparison with known exact results is also discussed. 
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I. INTRODUCTION 



Interest in strongly correlated electron systems and superconductivity has moti- 
vated some attention on the wide class of Extended Hubbard models |I[ . For such 
models, a certain number of exact results has been obtained for specific values of the 
interaction parameters 0,0, expecially in the one-dimensional case ||. These 
results support in particular the existence of a superconducting phase at T = 0, 
which is worth investigating also when relaxing some of the constraints among the 
parameters that allow the integrability of the Hamiltonian. 

The Extended Hubbard Hamiltonian reads 

Heh = H Hub + H x + H x + Hy + H w + H Y , (1) 

where 

H H ub = ~t E<y> E CT c U c j,a + u J2i n i,r n i,l (Hubbard) 

H X = X J2<ij> + n i,-«) c L c i,* (bond-charge) 

H x = X £ <i j > Yla n i,-o n i-o c \ a c j,<j (correl. hopping <-> on-site number) 

Hy = \ S<ij> n i n j (neighboring site charge) 

H w = t E<i,j> Ea,<T' c \A,°' c W c i,° (exchange) 

H Y = Y J2<^> c l,T c U c ja c j,T (P air hopping) . 

Here c\ a and c- x a are fermionic creation and annihilation operators, the sub- 
script i running over the sites of a <i-dimensional lattice A of L d sites, and 
° £ {T> 1} being the spin label; the usual anticommutation rules {c i a ,, Cj a } = 0, 
{ c i.o-> c \,a'} = ^ij ^c,cr' hold. The symbol < i, j > stands for nearest neighbors in 
A. Finally, n ia = c ia c ia . is the electron number operator with spin a at site i, and 
n, = 7i iT + 7i u . 

As Hubbard himself pointed out in [|J , all the terms in (|l|) appear due to the second 



quantization formulation of the interaction among electrons. The specific physical 
system of interest then suggests which are the most relevant ones. 

For the pure Hubbard model no exact result proves the existence of a super- 
conducting phase at finite values of U, and even within a mean-field scheme this is 
achieved only for negative U values. On the contrary a number of interesting exact 
results on the existence of a superconducting phase for appropriate non-vanishing 
values of some other interaction parameters in ([!]) has been obtained. Most of these 
results involve the notion of the states known as 7^-pairs, i. e. the states 

\v)4> = (Klm , m = l,...,L d (2) 

with 

jeA keB 
where is a <i-dimensional vector (<f), <ft, . . .), and B is the first Brillouin zone in the 
reciprocal lattice. Noticeably, the states [r])^ enjoy the property of 'Off Diagonal 



Long Range Order' (ODLRO), which implies superconductivity fL0 |. Therefore 
much effort has been done through the last years to find which are the relations 
among the coupling parameters in (JH) guaranteeing that an 1 77) ^ is the ground state. 
A first set of remarkable results in 1-d was obtained for the subclass of Hamiltonians 
characterized by the constraint X — t. In |6|,[j]] the phase diagram U vs filling for 
the AAS model X = t,X = W = V = Y = (reported in fig.0) was derived 
at T = : one can see a superconductive filling-independent region, where the 
l^)^ are degenerate ground states for any 0, and a filling dependent zone (again 
superconducting because it contains at least \r]) pairs) rising up to positive values 
of U. Unfortunately, in contrast to the real case of superconducting materials, 
the superconducting phase turns out to have a maximum at half filling. A similar 
phase diagram (see fig. 1) was also obtained in ]4|,|5| for the EKS Hamiltonian, 



characterized by X = t,X = 0, Y = W = V = — 1. There the filling independent 
phase is made of \rj)^ with just = 0, since Y ^ 0. In fact a non-vanishing pair- 
hopping term removes the degeneracy of \f])(p, \rj)o being energetically favourite for 

Y < 0, while \rj) v is favourite for Y > 0. Moreover, as fig.l shows, a non-vanishing 

Y also contributes to rise up the superconducting region towards positive values of 
U. 

More recently, it has also been realized || that, at least in order to obtain a 
filling independent superconducting ground state region, some of the constraints on 
the parameters (in particular X — t, which is not very physical) are not necessary, 
provided that a pair-hopping term is present (Y ^ 0). Notice that, as soon as 
I/i, only \r]) n states could become ground states, the other choices of in @ 
giving states that cannot be eigenstates of ([!]). 

On the contrary, no exact result is available concerning the existence of the more 
structured filling dependent superconducting region when the constraints on the 
parameters which allow the 1-d integrability are removed. It is one purpose of the 
present paper to investigate within the Hartree-Fock scheme such possibility, as well 
as to test how the superconducting region modifies at T ^ 0. 

Given the relevance of the pair-hopping term to stabilize the r^-pairs phase, and 
in order to make the physical mechanism more clear, we shall focus on a subcase 
of the extended model in which, apart from the pure Hubbard terms, only the pair 
hopping amplitude is taken different from zero. This is known in the literature as 
Penson-Kolb-Hubbard model. We want to emphasize here that the presence of other 
terms in (P is not expected to affect our results in a qualitative way, as other recent 



numerical studies confirm [p] ,p^|. 

In Section || we give the Hamiltonian and derive within the Hartree-Fock scheme 
the temperature dependent equations for the filling and the self-consistent supercon- 



ducting order parameter. In Section |IT| we solve the equations in dimension 2 and 
3, and show the temperature and filling dependence of the superconducting phase 
in these cases. In Section ^ we add to the Hamiltonian a next-nearest neighbors 
contribution to the hopping term, and we show how this affects the filling and tem- 
perature dependence of the superconducting phase. Finally, in Section [V] we discuss 
our results and give some conclusions. 

II. PENSON-KOLB-HUBBARD MODEL AND HARTREE-FOCK (HF) 
APPROACH TO THE 7T-PHASE 

The Penson-Kolb-Hubbard Hamiltonian reads 

hpkh = -t J2 c L c i,° + u J2 n ^ ni 'i + Y Yl 4t c u c ja c j,T + vJ2 n i > 

<ij> a i <i,j> i 

(4) 

where the last term accounts for the chemical potential. 

The case U = in l-d was first examined by Penson and Kolb |13[ to study a short 
range interaction between electron pairs of small radius (actually zero), in contrast 
with the BCS-theory, where the size of the Cooper pairs is comparatively large. 
This led to envisage a real space formulation for the electron pairing, which is very 
interesting in many contexts of condensed matter physics. 

Later on, the Coulomb repulsion term U was also taken into account by [14], where 



the PKH Hamiltonian was proposed as an effective phenomenological model captur- 
ing the main physical features of doped materials, such as high T c superconductors. 
Indeed if we assume that, due to some (yet unknown) microscopic mechanism, local- 
ized pairs can be formed, then their displacement in the lattice should be described 
by a pair hopping term competiting with a single carrier hopping amplitude t. The 
Coulomb repulsion term should account for the insulating phase. 



More recently a slave boson wide study of the different possible phases of the model 



at zero temperature has also been done |T5|]. In particular, a region characterized 
by a non-vanishing value of the order parameter = < K\ > was found; here 
K\ is given by (^) with <p = tt, and < > stands for the average value on the grand 
canonical statistical ensemble. In the following we shall denote such a phase as 
7r-phase. 

Noticing that H Y in (|l|) can be rewritten also as — Y X!<i j> ^Mt c n c j,l c j,t' 
and assuming a translational invariance throughout the lattice, a Hartree-Fock lin- 
earization can be implemented on both interaction terms in (fj) to study the 7r-phase. 
We shall suppose that such a phase is the energetically favourite one; therefore we 
shall take Y > and assume that the thermal energy is not sufficient to let other 
phases emerge (see |j5j for a detailed study of the possible several phases). Within 



this HF approximation the PKH Hamiltonian decomposes into a sum of /c-space 
independent hamiltonians, 



tt rp 

Hpkh ~ Y< H * ' [#k,#k']=0 , (5) 



HF 

where 

Hk = - [(ie k + ^)n kT + (te^-k + ju)n7r_kj.] + U ^7rC k T ct 7r-k| + ^ c 7r-kiC kT ^ - Ulx^ 2 

(6) 

and U = U — qY, q being the number of nearest neighbors, equal to 2d for a 
hypercubic lattice. The sum in runs over the Brillouin zone B and the k vectors 
are measured in units of the inverse lattice spacing (i.e. — tc < hi < 7r). 
In contrast with (§), the linearized hamiltonian (g) does not preserve the number of 
particles; indeed in a HF picture the 7r-phase has to be thought of as a superposition 
of ^-pairs involving different number of pairs, the average number of electrons being 
fixed through the chemical potential //. 



Since the dynamical algebra of each H k is the 577(2) © £7(1) given by the following 
generators: 



J f - c f c 1 " 

J k — C kT C 7T-k| 



7 k — c-7r_ k jC k | 



7 k = (riTT-ki + n-TT-kl - l)/2 
= ( n 7r-kT - «7r-k|)/2 
we can diagonalize it by means of the following unitary transformation, 



H' = e L *H k e- L 



, Lk — Ok( a; 7r^k ~~ X n^k 



where ak are real numbers satisfying the relation 

\Xtt\U 



tan (2ak|a; 7r | 



M + *(ek + e7r-k)/2 

The grand partition function Z in the HF approximation can then be worked out, 
since 



Z = Tr e ~ mpKH 



\{Tt^ h ^\{Tt^ h - 



where 7>k is the trace over the space of quantum labels (k, |) and {tz — k, {). 
The thermodynamic grand potential (per particle) in the thermodynamic limit now 
reads 



u = lim —ksT ■ „ 

T.D. N 



(7) 



(27T) 



rfk 



with 



A k + A; B T ln[2 cosh p D ** R * ] + ln[2 cosh /? Dk - ^ k ] 



S'k = t(e k + e 7r _ k )/2 
£>k = t(e k - e7r_ k )/2 
A k = /i + .Sk + t/|x^| 2 



(8) 



In order to investigate the thermodynamical properties of the system one has to 
implement the selfconsistency equation du)/dx n = for the order parameter, which 
turns out to be 



dk 



B 



1 + -7^- [tanh (3 — k ^~ ^ k — tanh/?- k k 







4R k ' 2 2 
As usual, this equation gives a solution x n = for T > T c and a solution x^ ^ for 
T < T c . It is easy to see that the 7r-phase exists only when U = U — qY < 0, that 
is when the pair-hopping Y > term renormalizes the effective interaction U > 
to an attractive regime. Indeed for T <T C the above equation gives: 



u- 1 



(2k) 



dk- 



B 



[tanh/3 — k -^ — - — tanh (3- k 



R\, 







2 '2 
Moreover, one must also satisfy the filling equation, n = du/dfj,, which reads 
1 



(9) 



n = 1 + 



(27T) 



dk [tanh ^^±^1 _ tanh (5^— ^] 

2i? k L P 2 P 2 J 



(10) 



Eqns. ©-([TOD constitute the parametric form of the equation of state. In order to 
get to one closed form, one should in principle invert (|T0"D obtaining as a function 
of n, T and {7, and then insert it into (0). The thermodynamics of the model will 
then be expressed in terms of n, T and U. 

Noticeably, comparing @ and (fTCf ) it is easy to show that whenever Sk in @ 
vanishes, z.e. whenever: 

e-7r-k = -ek (11) 
the chemical potential can be exactly inverted. In this case we simply have: 

»= l -^U (12) 

We wish to stress here that for a given model, even within the HF approximation, 
it is not obvious at all that /i can be inverted exactly: in the PKH model it is a 



peculiar feature of the 7r-phase (which is not shared by |r/)o phase). 

The relations (|i~l"l)-(|T^) hold in any dimension for a hypercubic lattice when dealing 

with a nearest neighbors hopping term, since the dispersion relations reads: 

d 

ek = 2 cos ki (13) 



As we will discuss in [IV|, eqn. (|TT|) is no more satisfied when a next nearest neighbors 
term is included. 



III. THE TEMPERATURE AND FILLING DEPENDENCE OF THE 
PHASE DIAGRAM IN D = 2 AND D = 3 

In this section we consider the case of a nearest neighbors hopping term, the dis- 
persion relation being given by (|HJ). We aim to derive the features of the critical 
values of U versus n for a given temperature T, in order to compare them with 
known solution of similar models (see section [V|). Such a critical curve U c = U c (n) 
is obtained from (|9|)- ([T0|) by setting the order parameter x n = into (H). 
As already noticed in section |I|, besides the parametric form (P)-(|10|), in this case we 
can also deal with one closed form; indeed, since ((TTJ) is satisfied, is independent 
of k, and therefore ( |12"D holds. By substituting /i into R^, and then R^ into (|), we 
obtain from ( |i0|) the critical equation 

1 



dk 



tanh Ut + U^n)^ _ ^ ^-^f(n)/2 



5(n) (14) 



(27r) rf y B 2 

where 5(n) = 1 — n is the 'doping' and B is the ci-dimensional Brillouin zone. 

In order to work out the integration it is worth to introduce the 'density of states' 

g^ d \t) with normalization 

WrL^L dt9i " {e) ' (15) 



which turns out to be in this case an even function of e, so that e m j n = — e max (this 
will not be the case in section |V|). Moreover we have e max = 2d. 



In d = 2 it is known [16] that 

^s^Hi)') ' (16) 

where K is the complete elliptic integral of the first kind; this function has a loga- 
rithmic divergence at e = 0. 

For the three-dimensional case we have also performed the calculation, obtaining 

2 f n 

- g (2) (e/2 -cos(x))dx < |e| < 2 

Jo 

g®(e) : = (17) 

2 />arccos(^— 2) 

-/ " g (2) (e/2 -cos(x))dx 2 < |e| < 6 , 
7T Jo 

which is plotted in fig.^|. 

Since eqn. (|14D is invariant under the transformation 5(n) — > — 5(n), it is easily seen 
that [/ c is symmetric with respect to the value at half filling (n — 1). The critical 
curves in d = 2 and d = 3 are plotted in figs.|] and |] respectively. For the sake of 
consistency with our approximation, we have plotted the region of the phase dia- 
gram where the values of U do not exceed the band- width A = 4td. 
Notice the different behaviour of U c in the two cases, in particular for low tempera- 
tures. Indeed in d = 2 we have a very sharp, cuspid-like shape at half filling, while 
in d = 3 a 'plateau' is obtained, meaning that the effective interaction threshold is 
almost independent of the density of electrons in the lattice for a rather wide range 
of n. This effect is due to the quite different behaviour of the density of states 
and g( 3 \ Indeed using the parametric form it is possible to show that at low T's the 
shape of U c in the neighborhood of half filling is governed by the behaviour of g^ 



around e = 0. In fact from eqns. (0) and (14) we easily obtain that the value of fi 
is determined by the equation 



gM(e)Pp{e;ri = l 

'^■m.n.x 



11 



where (3 = l/k B T, and Pp{e; /i) is 

Pp(e; //) = ^[tanh/%^-^ - tanh/%^-^] . 

Taking into account the normalization ([Tj^) of g^ d \e) and the fact that P^(e;/i) is 
positive and not greater than 1, one can easily see that n = 1 implies fi = (for any 
temperature T). Since from eqn. (|T2"D we have U = 2/i/(l — n), a simple limit — > 
can be performed, yielding 

-i 



U c (n = 1) = - 



<? (d) (^(e) 



(19) 



where 



/> ,(. i = — cosh 2 (^e) 



/ft 

T 



(20) 



This holds for any temperature T. In particular when T — > one gets 

1 



C/ C (n = 1) = - 



<?W(0) 



(21) 



Now, since g( 3 ) is almost constant around e = 0, when T ~ eqn. ( |T8| ) gives 



n 



g 



(3), 



le=0 



(22) 



and using ([TJ) we obtain that U c is actually independent of n (indeed g^(0) = 
0.1427 and so 1/0.1427 = 7.0078 which is just the value of U around half filling 
when T = 0). To be more precise, since the value of g( 3 ' at |e| = 2 is slightly higher 
than at e = (see fig.0), we have that at T = the plateau is slightly concave, 
so that the highest value of U c is actually reached away from half filling, at the 
symmetric values n ~ 0.4 and n ~ 1.6. 



Eqn. (|21|) stems from eqn.(|T^) because the functions pp(e) are 'regularizing' in the 



sense that they are positive and their integral on 1Z equals 1 for any (3 (though 
strictly speaking they belong to S{1Z) and not to C^ilZ)). In fact in d = 2 we 
observe that, when (3 — > oo, 

4 ,2/(00 

> / Pp(e) ■ inf g® oo , (23) 

7-2/09*) 0<|e]<2/C9t) 

V v ' 

2 tanh 1 

because of the divergence of g( 2 ' at e = 0. 

The 3-d case is slightly different because g^ is finite at e = 0; nevertheless, if we 
take a sufficiently small 5 around e = 0, it easily realized that 

Jim / g ( - 3) (e)p p (e)de = = Jim / P(s(e)de , 

J 5<\e\<e max ^°°Js<\e\<e max 

which yields 

g {3) (e) pfl(e) < sup # (3) V5 > , 



and so again (|21|) 



IV. THE NEXT NEAREST NEIGHBORS CONTRIBUTION 

Let us now turn to the case when in the Hamiltonian (f|) a next nearest neighbors 
contribution is included in the hopping term, which therefore becomes 

^EE c U,- at E E c U,* > (24) 

<i,j> cr <<i,j>> a 

Notice that such a term explicitly breaks the particle-hole cj — > e l7V '-*c^ invariance of 
the model. The dispersion relation is now given by 

d 

tk = 2 cos ki + a Y[ 4 cos k { cos A;^ . (25) 

i=l i<j<d 



As observed at the end of section |]J the symmetry = — e^-it does not hold 
anymore; as a consequence, Sk in eqn.(||) is no more vanishing, and does depend 
on k. This yields both mathematical and physical new features. In particular, the 
equation for the critical surface has to be given in the parametric form ©-(0), fi 
being the parameter. 

We are here interested in the case d = 2. We have computed also in this case the 
density of states, that turns out to be 



where a is assumed to be |a| < 1/2, since the next-nearest neighbors term is ex- 
pected to be small (in L 2 ) with respect to the nearest neighbors one. For a = we 
recover the usual form (^). Notice also that the function is no more even in 

("2") (2) 

e; indeed we now have ga (— e) = 9-a{ € )- The logarithmic divergence is shifted to 
e = —4a, and the domain is characterized by e m j„ = —4(1 — a) and e max = 4(1 + a), 
as shown in fig|| 

This imparts the critical curve U c vs n an asymmetric form, the highest value of U 
falling now at a n max ^ 1, as shown in fig.||. Moreover, at a given temperature, such 
a maximum of the critical curve is shifted upward with respect to the curve of the 
case a = 0. This means that at a given temperature T, the effect of the next nearest 
neighbors term is to reduce the 'optimal' effective attraction U. In turn, this implies 
at a given U the raising of the highest critical temperature reachable by doping the 
system. 

We have also studied how n max depends on a. The relation is almost linear for 
—0.4 < a < 0.4, while it displays a sudden increase around a ~ 0.45; in fig.0 we 
have extended the curve to the range < |a| < 1 (which could be still acceptable in 
principle) to show how n max saturates towards the limiting values n = or n = 2. 
Notice that the curve is odd in a; this is because the parametric eqns.fl9D-(|i~0D are 




(26) 



invariant under the transformation /i — > — /i ; a — > —a ; U — > [/, and therefore one 



can show, in agreement with [16|, that 



T c (U, n; a) = T c ([/, 2 - n; -a) . 

It is possible to see that the curve depends very weakly on the temperature T. 
Finally, we have plotted the phase diagram of T c versus n at fixed a, which is 
shown in fig.[8[ We can observe that the next nearest neighbors term has mainly two 
effects. First it shifts away from half filling the range of values of n at which the 
superconducting phase exists; this suggests that the system has to be doped in order 
to observe a superconducting behaviour. Secondly, it raises up the highest reach- 
able critical temperature with respect to the case where only a nearest neighbors 
interaction is considered. Having in mind the phenomenology of high T c materi- 
als, this study supports the idea that the actual microscopic Hamiltonian should be 
particle- hole not-invariant. 



V. DISCUSSION AND CONCLUSIONS 

The phase diagrams obtained within the Hartree-Fock scheme can be given a 
more precise physical interpretation by comparing them with the 1-d exact results 
known for some very specific cases. Strictly speaking, such a comparison is only 
possible for T = 0, since the phase diagram of these integrable models is not known 
at T ^ 0. In so doing, we observe that at zero temperature our results in fig.§ 
and fig.|] have the same structure as those in fig.|l[ One can recognize three distinct 
regions in the phase space. The first one is characterized by filling independent 
7r-phase. By reducing the effective attraction U, one enters a second region in which 
the existence of the 7r-phase depends on the actual filling. Finally, above the critical 
curve the 7r-phase disappears. 



In analogy with what is known to happen for the solved 1-d cases, we can think 
of the filling independent region as the phase in which all particles are paired in 
r^-pairs (Q) . This is in agreement with known result in d > 1 . The second region 
should be characterized by simultaneous presence of paired and unpaired electrons 
and empty sites, whereas in the third one no paired electrons could move. 

Switching on the temperature, thermal fluctuations are expected to break pairs. 
The dependence on the temperature of our phase diagrams supports this idea. Ac- 
tually for a given filling n, the greater becomes the temperature T, the greater must 
be the magnitude of the effective attractive interaction U in order to keep the re- 
pairs bound together. In fact this is what our phase diagrams in figs. [3] and [| show: 
the curves of higher T's lay below the lower T's ones; this can also be demonstrated 
analitically very easily, at least in the cases of a nearest neighbors hopping term. 
Indeed from eqn. fllif) , which is of the form G(n, (3, U c ) = 1 — n, we obtain that 
dUc l_d£c 1 dG/d(5 

dT k B T* df3 k B T*dG/dU c ' 1 J 

By working out the derivatives, and using the fact that the functions tanh(x) and 
xj cosh 2 (x) are monotone and that U is negative, it is easily seen that (p7|) is al- 
ways < 0. 

It is important to recall that, thanks to the presence of the pair-hopping term, an 
effective attractive interaction U is consistent with a positive value of the Coulomb 
interaction U. Hence the present HF treatment of the thermodynamics of the PKH 
model yields a structured filling dependent superconducting phase even in presence 
of repulsive on site Coulomb interaction between electrons. 

Finally, we wish to stress that the phase diagram in fig.| -obtained by including 
the next nearest neighbors hopping term- exhibits appealing features: for q^O the 
optimal doping of the superconducting region is at n max ^ 1, and the critical tem- 
perature is enhanced. Moreover, with respect to the results reported in jrjj on the 



attractive Hubbard model, our figure shows that even at T = the superconducting 
phase exists only for an appropriate range of filling values, not including half-filling. 
In fact, one could observe that the whole curve of the critical temperature vs. filling 
(HD is very reminiscent of that obtained for high-T c materials. 

The study of the influence that particle-hole non symmetric terms in the Hamil- 
tonian have on the features of the phase diagram has been worked out in 2-d, the 
conduction in high-T c superconducting materials tipically taking place along the 
cuprate planes. As Hartree-Fock approach is more accurate the higher is the dimen- 
sion, dealing with a 3-d and anisotropic order parameter would possibly be more 
reliable. Work is in progress along these lines. At the same time, since the results 
obtained here are encouraging, a numerical study of the temperature behavior of 
the present model in d — 2 would be probative. 
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FIGURES 

FIG. 1. Phase diagram in d = 1 at T = of AAS model (solid curve) and EKS model 
(dot-dashed curve). A filling-independent (superconducting) region is obtained in which 
the rj a pairs are ground state. Above it, a filling-dependent region (still superconducting) 
exists. In the EKS model the presence of a pair hopping term contributes to extend the 
superconducting zone towards positive values of U. In such exactly solved 1-d models 
the phase diagram reaches a maximum around half filling; moreover some unphysical 
constraint (for instance X = t) have to be imposed on the coupling constants in order to 
make the Hamiltonian integrable. 

FIG. 2. The density of states in d = 3 of the hopping term; notice that the value of 
g( 3 ) in the cuspids at |e| = 2 is slightly higher than at e = 0; this yields (see fig.f|) that at 
T = the maximum of the critical curve is reached away from half filling. 

FIG. 3. Phase diagram of the -zr-phase in d = 2; the critical value U c of the effective 
attraction U = U — qY is plotted as a function of the filling for some values of the 
temperature T. At a given temperature the -zr-phase exists for U < U c ; here A = 8t is 
the bandwidth and q = 4 is the number of nearest neighbors in the lattice. The values of 
U are negative because the pair hopping term renormalizes the on-site Coulomb repulsion 
U to a negative regime. Notice also that the curves are centered around half filling; at 
T = a filling independent region exists for U < —A, like in AAS and EKS models (see 



FIG. 4. Phase diagram of the 7r-phase in d = 3 ; A = 12t is the bandwidth and q = 6 is 
the number of nearest neighbors in the lattice. With respect to the bidimensional case (see 
fig.||) the curves have a plateau around half filling; indeed at T = the maximum values 
of U c are reached at the symmetric values n ~ 0.4 and n ~ 1.6. As the temperature is 
increased, the r/^-pairs become more likely to be broken and the extension of the 7r-phase 
reduces. 

FIG. 5. The density of states in d = 2 of the hopping term when a next nearest 
neighbors term of coupling constant a is also taken into account (solid line); notice the 
asymmetric behaviour with respect to the usual case a = 0: the logarithmic divergence is 
shifted to e = —4a. This effect is due to the break up of the particle-hole invariance of 
the hopping term; it yields an asymmetric shape of the phase diagram (see fig.|6|). 

FIG. 6. The 7r-phase diagram in d = 2 for different values of the next-nearest neighbors 
hopping amplitude a at the temperature ksT/t = 0.1. Notice that when a is increased the 
superconducting region rises up towards less negative values of the U = U — qY, and its 
maximum is reached at a n max which moves away from half filling. This means that, at a 
given temperature, the next-nearest neighbors interaction reduces the effective attraction 
U, yielding an increase of the highest reachable critical temperature (see also fig.||). 

FIG. 7. The behaviour of n max (i.e. the filling at which the critical curve U c vs n in fig|] 
reaches the maximum) as a function of the next nearest neighbors coupling constant a (in 
units oft), at temperature ksT/t = 0.1. The curve is odd. Notice that for —0.4 < a < 0.4 
the behaviour is almost linear; cis \cx\ is furtlierly incr63>s6d., Timax saturates to or to 2, as 
one can also see in fig.||. It is possible to show that this behaviour of n max depends very 
weakly on the temperature T. 



FIG. 8. The critical temperature as a function of the filling in d = 2, for a given value 
U = —4t of the effective attractive interaction. Here a is the coupling constant of the next 
nearest neighbors hopping term, which breaks the particle-hole symmetry of the model. 
Notice that, with respect to the case a = 0, the next nearest neighbors term yields both 
the increase of the highest critical temperature and the displacement away from half filling 
of the superconducting 7r-phase. In fact the highest T c is reached at n ~ 1.3. 
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